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We present a theoretical description of the diflerential conductance of point contacts between a 
normal metal and a multiband superconductor with extended pairing symmetry. We demonstrate 
that the interband impurity scattering broadens the coherent peak near the superconducting gap 
and significantly reduces its height even at relatively low scattering rates. This broadening is 
consistent with a number of recent experiments performed for both tunnel junctions and larger 
diffusive contacts. Our theory helps to better evaluate the energy gap of iron-based superconductors 
from point-contact Andreev spectroscopy measurements. 
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Introduction. Iron-pnictide superconductors are 
multiband materials with several disconnected Fermi 
surfaces. [ll-Q The structure of the superconducting (SC) 
energy gap near the Fermi surface (FS) remains one of the 
unsettled questions of the current investigation of these 
materials. Experiments indicate that in some compounds 
the SC gap vanishes along nodal lines, while many 
other observations are consistent with a fully gapped FS 
in SC states. 0-Q But even fully gapped SC states may 
be realized by different possible structures of the SC or- 
der parameter (OP), such as conventional or extended 
symmetry pairings, [li [loj 

Various experimental techniques have been utilized 
to explore the structure of the SC OP in iron 
pnictidcs, includin g a ngle-resolved photoemission spec- 
troscopy (ARPES)^, NMrEI and the London pene- 
tration depth technique. l^, 14 1 More recently, 
several groups have reported on measurements of the dif- 
ferential conductance of a point contact between a single- 
crystal iron-based pnictide superconductor and a normal 
metal by a technique called point-contact Andreev spec- 
troscopy (PCAS). Point contacts with both high tun- 
nel barrier s 11 5l4l7l and a large number of conduction 
channels p. [tI. Il8l42d | were reported. In either system the 
differential conductance curves were significantly broad- 
ened when compared to curves for conventional supercon- 
ductors within the Blonder- Tinkham-Klapwijk (BTK) 
theory pH and required the introduction of a somewhat 
short phenomenological electron lifetime (itI. [2^ near the 
FSs. Previous theoretical analysis pjj of PCAS for clean 
multiband superconductors does not seem to describe the 
above experiments. 

In the present Rapid Communication we argue that the 
smearing of the differential conductance in PCAS may 
indicate the realization of an extended s^-wave pairing 
and moderate interband impurity scattering. For super- 
conductors with extended s^-wave pairing, interband im- 
purity scattering mixes quasiparticle states between FSs 
characterized by opposite signs of the OP and has a pair- 
breaking effect on superconductivity. 13, 2^ A relatively 



small rate r^r ^ Tc of the interband scattering (we use 
h ^ 1 and fcs = 1) does not significantly change the crit- 
ical temperature Tc or the SC OPs A^.h on electron and 
hole pockets of s^-wave superconductors. We show, how- 
ever, that even weak interband scattering significantly 
modifies the excitation spectrum near the FSs and re- 
duces a resonance peak in the differential conductance at 
bias near the SC gap. 

This mechanism of broadening the differential conduc- 
tance peaks for s*-wave superconductors does not in- 
volve a phenomenological parameter for the electron life- 
time. The same model[25[ was applied to the explana- 
tion of the temperature dependence of the London pen- 
etration depth. 0, 0] The magnetic penetration depth at 
low temperatures exhibits a weak nonexponential depen- 
dence on temperature in electron-doped BaFe2As2 mate- 
rials, consistent with the theoretical temperature depen- 
dence for an extended SC state in the presence of in- 
terband impurity scattering. Interband impurity scatter- 
ing also results in a power-law dependence of the spin re- 
laxation rate in nuclear magnetic field resonance. [3, [l^l 
Therefore, one can evaluate SC OPs and interband scat- 
tering rates from independent experiments. 

Model. We consider a simplified two-band model[23l of 
iron-based pnictide superconductors with a single elec- 
tron and hole FSs. Generalization of this model to a 
larger number of electron and hole FSs does not qual- 
itatively change our results. We assume that the SC 
state has an extended s^-wave symmetry and is char- 
acterized by isotropic OP on electron, a = e, or hole, 
a = h, Fermi surfaces, with A/jAe < 0. Scattering off 
disorder can be separated into two categories: intraband 
scattering when the electron band index stays the same 
and interband scattering when the electrons change their 
band index. We note that the intraband scattering docs 
not affect the SC state with isotropic OP. The interband 
scattering mixes electron states with opposite OPs in 
wave SC and results in the depairing of the Cooper pairs 
suppressing superconductivity. In this case, the inter- 
band scattering rate F^ characterizes SC properties in 



2 



a similar way to the spin-scattering rate in conventional 
superconductors with magnetic impurities. (27| 

For a weak enough current through the point con- 
tact, the superconducting state is nearly homogeneous 
in space. We disregard the proximity effect in a metallic 
probe 28 1 and apply the circuit theory [2^ to evaluate 
the differential conductance G{V) = dl/dV: 



Cai2Tu±eV) 
2 cosh^ u 



du. 



(1) 



Here 



Im^Q 



(Rea„)2 ((Re^„ + 1)^ + (Im^a)') 



(2) 



has a meaning of spectral current density and is written 
in terms of functions aa{e) and ^q(£) = [1 — (£)]/[! -t- 
a^(e)] of quasiparticle energy e. The function Za(x) is 
defined in terms of the transmission eigenvalues t^n^ of a 
point contact between states in normal tip and electronic 
states on Fermi surface a: 
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The functions aa{e) arc solutions of a system of the 
following two fourth-order algebraic equations: 

An 2n , o- 2r^(Qe -a/i)(l +a/,ae) , . 
Ae(l - ttg) -I- 2i£ae = 2 ' (^a) 

+ '^h 

A n 2n , 2r^(Qfa -ae)(l +0^0/,) 
Aft(l - a,J + 2i£ah = — ^ . (4b) 

A proper solution of Eq. Q is chosen from the condi- 
tion that behaves as iAQ,/[2(e + iFtt)] at jej |Aq| 
and Eqs. ([I])-© recover the normal conductance Gat = 

(eV27r) Ea.„*«"^ at bias eV > |Ae,/i|. We note that for 
Ftt = 0, Eq. (HI) describes the BCS-type superconductor 
and Eq. ([T]) recovers the BTK result. 

Results. We first analyze the differential conduc- 
tance for a superconductor with equal-in-magnitude OPs, 
|Ae,/j| = A. The left panel in Fig.[T]shows the differential 
conductance for a single-channel junction with very weak 
tunnel probability t'"-* ^ 1 between the normal tip and 
the superconductor. 

In this case the Andrecv reflection is suppressed as 
(i^"-*)^ and the differential conductance is proportional 
to the electron density of states (DOS) in a superconduc- 
tor, except at the very top of Andreev resonance peaks 
at bias eV = |Ae_/j|. As the ratio F^/A increases, we 
observe an evolution of the DOS from a BCS-type DOS 
with gap equal to A and a very high peak above the gap 
to a smaller gap and a reduced height of the peak. In 
particular, a relatively small F^r/A = 0.01 already dras- 
tically reduces the height of the peak. At F^/A ~ 0.5, 
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FIG. 1. (Color online) Zero-temperature normalized differ- 
ential conductance curves in the case = — Ae = A for 
different values of r^/A = 0.01, 0.1, 0.5, 0.7. Left panel: A 
single-channel tunnel junction with t'''^ = f'*^' = 0.01. Right 
panel: A diffusive junction with a representative sample of 200 
channels with transmission eigenvalues distributed according 
to Eq. (O with a mean (fi''') = (fi"') = 0.84. 



the gap completely disappears, as expected. [2 Tj We also 
note that the position of the maximum of the differen- 
tial conductance moves slowly to higher bias as F^r/A 
increases. 

Larger contacts have many conduction channels with 
transmission eigenvalues ti"^ between and 1. Assuming 
that the contact is diffusive, we model the statistics of 
in"'' by the Dorokhov distribution: jsof 



P(e)) 



1 



1 — 



(a) 



sech^Qc 



<e^<l (5) 



with a relatively high cutoff qa 1. The parameter 
determines the average (in"'') = tanhgti/ga, and the total 
number of channels can be chosen to match the junction 
conductance in the normal state Gat. We present the 
differential conductance for such diffusive contacts in case 
A/j = — Ae = A for several ratios of F^/A in the right 
panel of Fig. [T] Unlike G{y) of tunnel junctions, G{V) of 
diffusive contacts may exceed On due to the contribution 
to the current from the Andreev reflection that doubles 
the conductance of nearly open channels with t^'^ ~ 1. 
As interband scattering rate F^ increases and suppresses 
superconductivity, G(y) curves smoothen and approach 
Gn- In the gapless regime, F^/A > 0.5, the differential 
conductance exhibits a zero-bias peak. 

In Fig. [2] we present a fitting of recent experimental 
data for diffusive point contacts using our model of 
an superconductor with equal OPs on the two FSs. 
Overall, large conductance indicates that there are many 
channels across the contact between the metallic tip and 
the (Bao.6Ko.4)Fe2As2 superconductor (Tc ~ 37 K). For 
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FIG. 2. (Color online) A fit to the experimental data of 
Ref. [ii, Fig. 3], with Ah = ~Ae = A = 6.3 meV and 
= 0.59 meV used as fitting parameters. Only experimental 
points for V > were utilized for the fit. The distribution of 
transmission eigenvalues, Eq. ((5|, was adjusted to match the 
experimental zero-bias conductance G(0). The temperature 
was fixed at the reported value T = 2 K. 




FIG. 3. (Color online) Zero-temperature normalized differen- 
tial conductance curves in the case F^/A = O.Of for differ- 
ent values of the ratio |Ae/Ah| = f.O, f.2, f.6 (A^ = A), 
(a) A single-channel tunnel junction with the sign of the OP 
opposite on different FSs. (Inset) A diffusive junction with 
a representative sample of 200 channels with transmission 
eigenvalues distributed according to Eq. ((5]), and with a mean 

(in'- ) = (tn ') = 0.84. The sign of the OP is opposite on dif- 
ferent FSs. (b) A single-channel tunnel junction with the sign 
of the OP equal on different FSs. 



an s superconductor with equal OPs |Ae.h| = A, the 
zero-bias conductance G{0)/Gn is independent from A 
and and depends only on the average value of trans- 



mission eigenvalues, chosen as (in ) = {tn ) ~ 0.84 to 
match the data. Then we found a good fit to experimen- 
tal data by taking the SC OP A = |Ae,h| = 6.3 mcV 
and = 0.59 mcV (see Fig. [5]). The OP A is somewhat 
greater than that found in Ref. [l^ based on a model 
with a phcnomcnological electron lifetime. (2^ We esti- 
mate 2A(0)/Tc « 3.95, which is close to the BCS value 
3.53 and consistent with enhancement of 2A/Tr. due to 
intcrband scattering in superconductors. [25j For the 
above parameters, we find r^/27rrc,o — 0.023, and the 
critical temperature is only moderately suppressed as 
Tc = 0.77 Tc^o due to the interband scattering when com- 
pared to idealistic hi the absence of scattering. We 
would like to emphasize that the value of F^ can be es- 
timated for the same sample from independent measure- 
ments of magnetic penetration depth or spin relaxation 
time. 

We also consider the differential conductance for an 
state with different magnitudes of OPs on the FSs, 
Ae/A/i ^ —1. For tunnel junctions with t*^"' <C 1, the 
differential conductance acquires a double-peak feature 
once the OPs on the two FSs are different [see Fig.[3l^a)]. 
We note that once the SC energy gaps on the FSs have 
different magnitudes, the resonance peak at smaller bias 
becomes sharper than that in the case of equal OPs. The 
sharp first peak is a consequence of suppression of the 
intcrband scattering of electrons, when one FS is still 
gapped while the other FS has quasiparticles. This be- 



havior is consistent with the observed differential conduc- 
tance of tunnel contacts in Ref. jT6| . As a comparison, 
Fig. El^b) presents the case of an s++ state characterized 
by equal signs of the OPs on the two FSs. Although the 
smea ring of the peaks is present even in this case (see also 
Ref. [3l| for DOS in an s++ state), the peaks are much 
sharper when compared to the state for |Ae/A;i| close 
to unity. When this ratio is away from unity, the G(V) 
curves for the and s^'^ states become similar, cf. plots 
for |Ae/A,,| = 1.2 and |Ae/A,,| = 1.6. 

In the inset of Fig. EJa) we present the differential con- 
ductance of a diffusive point contact. We note that the 
low-bias coherence peak does not change much when the 
OPs are taken with different magnitudes, since this peak 
is already broadened by a simultaneous contribution of 
the tunneling current above the gap and Andreev reflec- 
tion within the gap. But the double-peak feature of the 
plot still remains, and the overall differential conductance 
peak appears to be broader than that for a model with 
equal magnitude gaps. 

To investigate temperature dependence of G{V), we 
consider Ah = — Ae = A(r). In this case the self- 
consistency equation for the OP reduces to 
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FIG. 4. (Color online) Temperature dependence (the varying 
parameter is T/Tc — 0, 0.1, 0.5, 0.9) of the normalized differ- 
ential conductance curves in the case = — Ae = A(T) for 
a choice of the parameters such that rTT/2nTc = 0.029. Left 
panel: A single-channel tunnel junction. Right panel: A dif- 
fusive junction with a representative sample of 200 channels 
with transmission eigenvalues distributed according to Eq. ([5]) 
and with a mean (t'n^) = =0.84. 



with an auxiliary variable ym that satisfies 



(6b) 



and ^'(a;) being the digamma function. The transi- 
tion temperature Tc is suppressed by the impurity in- 
terband scattering with respect to the clean critical 
temperature Tc^ according to \n{Tc/Tcfi) = ^"(1/2) — 
^(1/2 + T^/ttTc), and Tc vanishes completely at T^r > 
0.07(27rTc,o). At T^ ~ 0.064(27rTc,o) a SC state be- 
comes gapless and has the critical temperature Tc ~ 
0.22Tc,o-U^ 

In Fig.|4]we plot G{V) for a tunnel junction and a dif- 
fusive point contact using Tc = 37K at T^ = 0.59 meV, 
found from the fit to experimental data in [loj with 
A(0) = 6.3 meV. We observe that in tunnel junctions, 
the positions of the conductance peaks are nearly in- 
dependent of temperature, but their heights decrease. 
The value of the zero-bias conductance grows monoton- 
ically with temperature until the curve flattens out at 
the critical temperature. In diffusive contact, the coher- 
ence peaks are rounded relatively quickly and move to 
zero bias as temperature increases. At T ~ 0.5 Tc the 
peak settles at y = 0, and further increase in tempera- 
ture leads to the suppression of this zero-bias peak near 
T = Tc. 

In conclusion^ we presented a theoretical description 
of the differential conductance of a point contact be- 
tween an iron-based superconductor and a metallic tip. 
We demonstrated that interband impurity scattering of 
quasiparticles between Fermi surfaces significantly mod- 
ifies differential conductance and provides a microscopic 
dcpairing mechanism in an s^-wave superconducting 
state. This behavior stems from the sign-changing su- 



perconducting order parameter between Fermi surfaces 
and is an indirect indicator of superconductivity 
in iron-based superconductors. The broadening of the 
differential conductance curves was observed in large 
contacts. Measurements of the differential conduc- 

tance of tunnel contacts [l5, 16 1 also showed a significant 
suppression of the coherence peaks in conductance at bias 
near the superconducting gap. These experimental ob- 
servations required an incorporation of a non-vanishing 
value of the depairing rate, which has a simple micro- 
scopic origin in our model. 

Electronic band structure and disorder in real materi- 
als are likely more complicated than those analyzed here. 
However, we expect that the relations between thermo- 
dynamic and transport characteristics of iron-based su- 
perconductors do not crucially depend on details of band 
structure and disorder and that a single parameter, such 
as the interband scattering rate T^y, describes them. Such 
a parameter can be evaluated from various independent 
experiments, such as PCAS and temperature dependence 
of the magnetic penetration depth that was also fitted 
in terms of . j25j Since impurity scattering rates often 
depend on conditions of material preparation, a more 
conclusive analysis would require comparison of results 
of PCAS and magnetic penetration depth measurements 
performed on samples produced in similar conditions, as 
was the case in recent studies of electron doped BaFe2 As2 
films. [3, [H 
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